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Abstract 

We present a numerical kinetic scheme for an unsteady mixed 
pressurized and free surface model. This model has a source term 
depending on both the space variable and the unknown U of the 
system. Using the Finite Volume and Kinetic (FVK) framework, we 
propose an approximation of the source terms following the principle 
of interfacial upwind with a kinetic interpretation. Then, several 
numerical tests are presented. 

Keywords : Finite volume scheme, Kinetic scheme, conservative source 
terms, non-conservative source terms, friction 

1 Introduction 

In this paper, we study a way to upwind the source terms of a mixed flows 
model in non uniform closed water pipes in a one dimensional framework. 
In the case of free surface incompressible flows, the model is called FS- 
model and it is an extension of classical Saint- Venant models. When the 
pipe is full, we introduce the pressurized model, called P-model, which 
describes the evolution of a compressible inviscid flow and is close to gas 
dynamics equations in a nozzle. In order to cope with the transition be- 
tween a free surface and a pressurized model, we use a mixed model called 
PFS-model. It is based on balance laws and provides an hyperbolic system 
with source terms corresponding to the inclination of the pipe (seen as a 
topography term) , the section variation, the curvature and the friction. 

Several ways to compute the numerical approximation of conservation 
laws with source terms have already been investigated by many authors. 
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The main difficulty is to preserve numerically some properties satisfied by 
the continuous model: the invariant domain, the well-balanced property for 
instance. The Finite Volume methods are largely used since they present 
the remarkable property to be domain invariant (for instance, for Saint- 
Venant equations, to be water height conservative). Some Well-Balanced 
Finite Volume schemes, introduced by Greenberg et al [3], preserve steady 
states. All these methods are based on two principles: firstly, the con- 
servative quantities are cell-centered as usual finite volume schemes, and 
secondly the source terms are upwinded at the cell interfaces. 

In this paper, we consider a particular Finite Volume- Kinetic scheme 
built to compute the numerical solution of PFS-model. This scheme is 
based on the classical kinetic interpretation [5] of the system. 

The source terms appearing in the PFS-model are either conserva- 
tive, non-conservative or else. All source terms are upwinded at the cell 
interfaces: we use the definition of the DLM theory [3] to define the non- 
conservative products. The particular case of the friction term which is 
neither conservative nor non-conservative will be upwinded using the no- 
tion of dynamic slope. The source terms are taken into account in the 
numerical fluxes and are computed from the microscopic ones, obtained 
through the concept of potential barrier. 

The paper is organized as follows. In the second section, we describe the 
PFS-model |2| and focus on the source terms. The detailed description of 
the method used to deal with the transition points (when a change of state 
occurs) is not presented (see [5] for more details on this topic). We state 
some theoretical properties of the system. In the third section, we give 
the kinetic formulation of the PFS-model and the corresponding kinetic 
scheme. In the fourth and last section, several numerical tests are provided. 

Notations concerning geometrical quantities 

• 6(x): angle of the inclination of the main pipe axis z = Z(x) at position x 

• Q(x) = S(x): cross-section area of the pipe orthogonal to the axis z — Z(x) 

• R(x): radius of the cross-section S(x) orthogonal to the axis z = Z(x) 

• £l(t, x) : free surface cross-section area orthogonal to the axis z = Z(x) 

• a(x, z): width of the cross-section Q. at position x and altitude z 

Notations concerning the PFS-model 

• p(t, x,y, z): pressure 

• pa: density of the water at atmospheric pressure po 

• p(t,x,y, z): density of the water at the current pressure 




p(t,x, y, z) dy dz: mean value of p over Q(x) 
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• c: sonic speed 




S(x): equivalent wet area 



• u(t, x): velocity 

• Q(t, x) = A(t, x)u(t, x): discharge 

• E: state indicator equal to E = if the flow is free surface, E — 1 otherwise 

• S: the physical wet area equal to A if the state is free surface, S otherwise 

• 'H(S): the ^-coordinate of the water level equal to "H(S) = h(t,x) if the 
state is free surface, R(x) otherwise 

• p(x, A, E): mean pressure over 

• K 3 > 0: Strickler coefficient depending on the material 

• Pm(A): wet perimeter of A (length of the part of the channels section in 
contact with the water) 

• Rh(A) = — — — — ■: hydraulic radius 

J m (A) 

• Bold characters are used for vectors, except for S 

2 A model for unsteady water flows in pipes 

The PFS-model [5] is a mixed model of a pressurized (compressible) and 
free surface (incompressible) flow in a one dimensional rigid pipe with vari- 
able cross-section. The pressurized parts of the flow correspond to a full 
pipe whereas the section is not completely filled for the free surface flow. 
The Free Surface part of the model is derived by writing the 3D Eu- 
ler incompressible equations and by averaging over orthogonal sections to 
the privileged axis of the flow. In the same spirit, by writing the Eu- 
ler isentropic and compressible equations with the linearized pressure law 

pit, x, y, z) — po H — oipitt x, y, z) — po), we obtain a Saint- Venant like sys- 
cr 

pit, x) 

tern of equations in the "FS-equivalent" variable A(t,x) — ■ — S(x), 



Q(t,x) = A(t,x)u(t,x) which takes into account the compressible effects 
(for a detailed derivation, see [2]). 

In order to deal with the transition points (that is, when a change of 
state occurs), we introduce a state indicator variable E which is equal to 
1 if the state is pressurized and to if the state is free surface. We define 
the physical wet area by: 



Po 




E=l, 
E = 0. 
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The pressure law is given by a mixed "hydrostatic" (for the free surface 
part of the flow) and "acoustic" type (for the pressurized part of the flow) 
as follows: 

p(x,A,E) = c 2 (A- S) + gh{x, S)cos9 (1) 

where g is the gravity constant, c the sonic speed of the water (assumed to 
be constant) and 8 the inclination of the pipe. The term Ii is the classical 
hydrostatic pressure: 



f H(S) 

h(x,S) = / (H(S)-z)adz 

J-R 



A 



where <r(x, z) is the width of the cross-section, R = R(x) the radius of the 
cross-section and 'H(S) is the z-coordinate of the free surface over the main 
axis Z(x). 

The defined pressure {T]) is continuous throughout the transition points and 
we define the PFS-model by: 

f d t (A)+d x (Q) = 

d t (Q) + d x ( %+p{x,A,E)\ = -gAZ' + Pr{x,A,E) 

J (2) 
-G(x,A,E) v ' 

where z = Z(x) is the altitude of the main pipe axis. The terms Pr, G 
and K denote respectively the pressure source term, a curvature term and 
the friction: 

Pr(x,A,E) = c 2 - l^j S' + gI 2 (x,S)cos6, 
G(x,A,E) = gAZ(x ) S)=gA(H(S)-I 1 (x,S)/S)(cos6)', 
K(x,A,E) ' 



i^ s 2 ^(S) 4 / 3 

where we have used the notation /' to denote the derivative with respect 
to the space variable x of any function f(x). The term I2 is the hydrostatic 

nU(S) 

pressure source term defined by: l2(x, S) = / (^(S) — z)d x adz . The 

J-R 

term K s > is the Strickler coefficient depending on the material and 
i?/j(S) is the hydraulic radius. 

The System §2§ has the following properties: 

Theorem 2.1 

1. System ([2]) is strictly hyperbolic on {A(t, x) > 0} . 
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2. For smooth solutions, the mean velocity u = Q/A satisfies 

d t u + d x f y + c 2 ln(A/S) + gH{S) cos6 + gZ^j (3) 
= -gK(x,A,E)u\u\^0. 

3. The still water steady state, for u = 0, reads: 

c 2 \n(A/S) +gH(S) cos 8 + gZ = 0. (4) 

4. System ([5]) admits a mathematical entropy 

O 2 — 
£(A,Q,E) = + c 2 A\n(A/S) + c 2 S + gAZ(x,S) cos0 + gAZ 

which satisfies the entropy relation for smooth solutions 

d t S + d x {{£ +p{x,A,E))U) = -gAK{x,A,E)u 2 \u\ s^O. (5) 

In what follows, when no confusion is possible, the term K(x,A,E) will 
be noted simply K(x, A) for free surface states and K(x, S) for pressurized 
states. 



3 The Kinetic approach 

The kinetic formulation (|3.ip is a (non physical) microscopic description of 
the PFS-modcl provided by a given real function \ '■ K — > K satisfying the 
following properties: 

xM = x(-w) >0 , x(u)duj = 1, / w 2 x(w)rfw = 1. 
Jw. Jm. 

It permits to define the density of particles, by a so-called Gibbs equilibrium, 
M(t,x,0 = ffi^x f ^T/^V^ where b(t,x) = b{x,A(t,x),E{t,x)) 



if £ = 0, 
if E = 1. 
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3.1 The mathematical kinetic formulation 

The Gibbs equilibrium M. is related to the PFS-model by the classical 
macro-microscopic kinetic relations: 

A= f M(t t x,£)dt, (6) 

Q = J tM(t,x,£)dt, (7) 
+ Ab(x,A,E) 2 = [ Z 2 M(t,x,Od£. (8) 
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From the relations ©-(jS]), the non-linear PFS-model can be viewed as a 
single linear equation involving the non- linear quantity M: 

Theorem 3.1 (Kinetic Formulation of the PFS-model) (A, Q) is a 

strong solution of System ^ if and only if Ai satisfies the kinetic transport 
equation: 

d t M+^d x M-g<t>d^M = K,{t,x^) (9) 
for a collision term K,{t, x,£) which satisfies for (t, x) a.e. 

The source terms are defined as: 

<f>(x,W) =B(x,W) -d x W (10) 



with 



W= (z + J K(x,A)u\u\dx, S, cosd) 

'yix.A^cosQ — . ,,\ 
1, -^-^ , Z(x, A)) if £ = 



(11) 



and B = 

I (^ 'BEl 

A 

where l2(x,S) reads j(x,S)S' for some function 7 (depending on the ge- 
ometry of the pipe). 

it is time and space variable dependent. 
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3.2 The kinetic scheme 

Based on the above kinetic formulation (0), we construct easily a Finite 
Volume scheme where the source terms are upwinded by a generalized 
kinetic scheme with reflections [5]. 

To this end, let us consider a uniform mesh on K where cells are denoted 

for every i G Z by mi — {x^_i/2, 2^+1/2)1 with Xi = 1 - — l+1 / 2 anc i 
Ax = x i+ i/2 — x i+ i/ 2 the space-step. We consider a time discretization 
t n defined by t n+x = t n + At" with At n the time-step. We note U^ 1 = 

(Af, Q"), u™ = — l -, M™ the cell-centered approximation of U = (A, Q), u 
A i 

and Ai on the cell rm at time t™. 

If W is {^Z + J K(x, A)u\u\ dx, S, cos d^j , its piecewise constant rep- 
resentation is given by, W(t,x) = Wi(t)l mi (i) where W;(t) is defined as 
^i(i) = -r- f W(i, x) for instance. 

Denoting by W,; and W, + i the left and the right states of the cell interface 
£1+1/21 an d using the "straight lines" paths (see [I]) 

*(s,W„W 4+1 ) = sW i+1 + (1 - »)W,, s e [0, 1], 

we define the non-conservative product <p(t, 2^+1/2) by writing: 

[W](t)- / B(t,*(a,W i (t),W < +i(t)))da (12) 
Jo 

where [W] (t) := Wj+i(i) — W,(t), is the jump of W(i) across the dis- 
continuity localized at x = ££+1/2 • As the first component of B is 1, we 
recover the classical interfacial upwinding for the term Z (appearing e.g. 
in Saint- Venant equations) since it is a conservative term. 

Neglecting the collision kernel [5] and using the fact that <fi — on the 
cell rm (since [W] = 0), the kinetic transport equation (j9|) simply reads: 

aa; , (t, x, € [i n , t„+i) xnijxR (13) 

> f(tn,X,0 =M(t n ,X,g) 

and thus it may be discretized as follows: 

A/" 

- + e (x~ + , (0 - M+_ 4 (0) (14) 



where the contribution of the source term is included into the microscopic 
1/2 



numerical fluxes jM i±1 / 2 - This is the principle of interfacial source upwind. 
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Using the macro-microscopic relations (J5])-© and integrating Equation 
(fT4|) against £ and £ 2 , we obtain the Finite Volume scheme: 

A/™ 

u ? +1 = u r + ^(F^-Fti) (16) 

where the numerical fluxes are computed by : 



(16) 



Following [5] (or pQ), the microscopic fluxes are given by: 

positive transmission reflection 



m- +1/2 (o = i u>0} Mm +t u<0 , e - 2gA ^ +i/2<0} Mu-o 



1{«<o,^-2 9 a^ +1/2 >o}M? +1 (^-y^ 2 - 2gA<% +1/2 ), 

negative transmission 
negative transmission reflection 



x++i /2 (o = i{f<o}^s+i(o +i{c>o^+2 9 A C+1/2 <o}xr+i(-e) 

+ l«>0,£*+2 9 A* r+va >0}A<? + 2. 9 A^ +1/2 ) . 

s « ' 

positive transmission 

(17) 

The term A0" ±1 ^ 2 in ([T71) is the upwinded source term ([TO]) . It also plays 
the role of the potential bareer: the term £ 2 ± 2gA<fiV' +1 , 2 is the jump 
condition for a particle with a kinetic speed £ which is necessary to 

• be reflected: this means that the particle has not enough kinetic 
energy £ 2 /2 to overpass the potential barrier (reflection in (|T71) )). 

• overpass the potential barrier with a positive speed (positive trans- 
mission in (fT7|). 

• overpass the potential barrier with a negative speed (negative trans- 
mission in (|17|l )). 

Taking an approximation of the non-conservative product cj) (|12p. the po- 
tential barrier A(/>™ +1 / 2 has the following expression: 

- [W] (t n ) • B (t n ,V (i, Wi(t n ), W i+1 (t„)^ (18) 

1 

2vT 



1/2 - L VT J V-n) ■ yb n , * y- 
Next, with the simplest choice of the %-function x{ u ) — 7T7k^\-^/3 V3] ( w )> 



which allows to compute easily numerical fluxes, we have: 
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Theorem 3.2 

1. Under the CFL condition — — max ( \u?\ + Vic ) < 1, the kinetic 

Ax iez V ' / 
scheme ([T5]) - ([TT1) keeps A positive, i.e. A™ ^ if it is initially true. 

2. The kinetic scheme (|15p - (fl7|) allows to compute the drying and flood- 
ing area. 

4 Numerical results 

Let us recall that the zero water level corresponds to the main pipe axis. 
The piezometric head (or line) is defined by: 

{ c 2 [A — S) 
p = 2R H if the flow is pressurized 

gS 
p = h if the flow is free surface, 

where h is the water height. 

Comparison with the VFRoe scheme |2j. 

We compare the result obtained by the presented kinetic scheme with 
the upwinded VFRoe method [2]. 

The numerical experiment is performed in the case of an expanding 5 
m long closed circular pipe at altitude Zq = 1 m with slope (slope of 
the main pipe axis). The upstream diameter is 2 m and the downstream 
diameter is 2.2 m. The friction is not considered for the first test and is 
set to 0. The simulation starts with a still free surface steady state. The 
upstream boundary condition is a prescribed piezometric line (increasing 
linearly from 1 m to 3. 2m in 5s) while the downstream discharge is kept 
constant to 0m 3 /s. The other parameters are N — 100 (discretization 
points), CFL= 0.8 and the sound speed is 20 m/s. 

The result is in a good agreement and is represented on Fig. [T] 



Upwinding of the friction. 

It is well-known that cell-centered approximation of source terms leads 
to, generally, wrong results. We consider the kinetic scheme with the up- 
winded friction and the cell-centered one (i.e. we use W = (Z, S, cos 9) 
instead of (|TT|) and we add the cell-centered friction K™u™\u™\ to the right 
hand side of Equation (|15[) V We compare the schemes in a symmetrical 
flow. 

The numerical experiment is performed on a 100 m long closed pipe with 
constant section of diameter 2 m. The simulation starts from a "double 
dam break", as displayed on FiG.[2]and FiG.[3]at time t = 0. The upstream 
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and downstream condition are identical: the piezometric head increases 
linearly from 1 to 2.1 meters. We choose the same parameters as in the 
previous experiment. 

The results in Fig. [2H3] show that the scheme with the cell-centered 
friction, contrary to the upwinded one, does not preserve the symmetry of 
the flow. In particular, for K s = 100 (low friction) and at time t = 56.210 
(see FiG.[2]on top) we observe a small disymmetry, which evolves drastically 
at time t = 1.095 for K s = 10 (high friction) (see FiG.[3]on top). Despite the 
unavailability of experimental data, the kinetic scheme with the upwinded 
friction term, from a physical point of view, gives the expected result, 
namely, a symmetrical flow. 

5 Conclusion 

We have presented a global manner to upwind conservative and non con- 
servative source terms. To this end, we have used the definition of the 
non-conservative product of [1] which allows to recover the classical up- 
winding of conservative terms. Using the notion of dynamic slope, we have 
also upwinded the friction term given by the Manning-Strickler law (which 
is neither conservative nor non-conservative) in a FVK framework. The 
combination of all these quantities into a single one is an elegant and easy 
way to construct a kinetic scheme with reflections by introducing the po- 
tential bareer. Although kinetic schemes naturally deal with drying and 
flooding areas, the friction term is manually set to when such cells ap- 
pear. 
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Figure 1: Comparison between the kinetic scheme and the upwinded 
VFRoe scheme at x — 0.5 m. 
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Figure 2: Comparison of the cell-centered friction and upwinded friction 
for K s = 100. 
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Figure 3: Comparison of the cell-centered friction and upwinded friction 
for K s = 10. 



13 



